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Abstract This paper discusses the observed flat rotation curves of galaxies 
in the context of noncommutative geometry. The energy density of such a 
geometry is diffused throughout a region due to the uncertainty encoded in 
the coordinate commutator. This intrinsic property appears to be sufficient 
for producing stable circular orbits, as well as attractive gravity, without the 
need for dark matter. 
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1 Introduction 

The inability to account for stellar motions in the outer regions of galaxies has 
led to the hypothesis that galaxies and even clusters of galaxies are pervaded 
by dark matter 1 2. '5 . This hypothesis has been confirmed by observing the 
flatness of galactic rotation curves [H0EH7JIMH] " ^ e composition of such 
matter, if it exists, has remained unknown. 

The dark-matter problem originated in the measurement of the tangential 
velocity of stable circular orbits of hydrogen clouds in the outer regions 
of the halo. To explain the observed constant velocity, it is assumed that the 
decrease in the energy density is proportional to r~ 2 , where r is the distance 
from the center of the galaxy. 

A number of candidates for dark matter have been proposed. The most 
favored is the cold dark matter (CDM) paradigm [TUlIll) . Another strong 
possibility is the yl-CDM model that is related to the accelerated expansion 
of the Universe p~2l[T3j . For a summary of some alternative theories, such as 
scalar-tensor and brane- world models, see Rahaman, et al. P^IITMTB] . 

In this paper we study the dark-matter problem from a completely differ- 
ent perspective: an important outcome of string theory is the realization that 
coordinates may become noncommuting operators on a D-brane 18, 19 . The 
result is a discretization of spacetime due to the commutator [x M , x"]= i 8^, 
where is an antisymmetric matrix. Noncommutativity replaces point-like 
structures by smeared objects [50], suggesting the possibility of eliminating 
the divergences that normally appear in general relativity. The smearing ef- 
fect is accomplished by using a Gaussian distribution of minimal length y/a 
instead of the Dirac-delta function. More precisely, the energy density of 
the static and spherically symmetric smeared and particle-like gravitational 
source has the following form (21) : 

, M ( r 2 \ 

where the mass M is diffused throughout a region of linear dimension y/a due 
to the uncertainty. The noncommutativity is an intrinsic geometric property 
of the manifold, and not of its matter content, and can be taken into account 
by keeping the standard form of the Einstein tensor on the left-hand side of 
the field equations. The right-hand side is modified, however, by introducing 
a new energy-momentum tensor as a gravitational source. 

Taking the flat rotation curves as input, this model predicts both stable 
circular orbits and attractive gravity in a typical galaxy. 

2 The model 

In this paper the metric for a static spherically symmetric spacetime is taken 
as 

ds 2 = -e^dt 2 + e x ^dr 2 + r 2 {dB 2 + sin 2 0d<f> 2 ), (2) 

where the functions of the radial coordinate r, v(r) and A(r), are the metric 
potentials. 
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Now we consider the model with a maximally localized source of energy. 
Here the Einstein equations can be written as 



(3) 



(We assume that c = 1). The most general energy momentum tensor com- 
patible with static spherical symmetry is 



T„ = diag{-p, p r , p u p t ). 
For the metric @, the Einstein field equations are 
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(5) 

(6) 
(7) 



3 The solutions 

Using the observed flat rotation curves as a starting point, it is well known 
15,22] that this condition gives the solution 

e v = B r l , (8) 

where I is given by I = 2v 2 ^ and Bq is an integration constant. (For a deriva- 
tion, see Ref. [IB].) According to Matos, Guzman and Lopez [53], the ob- 
served rotation curve profile in the presumed dark matter dominated region 
is such that the rotational velocity becomes approximately constant with 

~ 300 km/s (~ 10~ 3 ) for a typical galaxy. So I = 0.000001, as shown by 
Nandi et al. |24j . and we likewise assume large distances measured in kpc 
from the galactic center. 

Using equation (fT]), equation §5§ yields 

e"* = 1 - (9) 



where 



and 



, 2M /3 r 2 \ 2M f r2/Aa r tj 

m ^=Vt 7 {2^) = Vtl '" (10) 



o 



3 



r 2 /4c 







is the lower incomplete gamma function [3T]. The classical Schwarzschild 
mass is recovered in the limit as r/y^a — > 00. Furthermore, new physics can 
only be expected if r ss y/a. 
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Fig. 1 The graphs of m*(r) (left) and m(r). 



The mass M could be a diffused centralized object such as a wormhole 
[2"5] or a gravastar [35]. Since we are interested in rotation curves at some 
fixed distance r = Rq from the center, we need to consider instead a spherical 
shell of radius r — Rq. So instead of a smeared particle-like object, we get a 
smeared shell. We will therefore need to translate the above curves as follows: 



M ( (r - R ) 2 

P = — TT ex P 



4o 



and 

< \ 2 

m\r\ = — , 

R V 2 4 « 

These replace p* and m*(r). As before 



2M /3 (r-i? ) 2 
m ( r ) = —7= 7 o> 



(12) 



(13) 
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(r-iio) /4a 



(14) 



is the lower incomplete gamma function. It is important to note that Eq. (|T4|) 
consists of a pure translation of Eq. (jllj) and therefore remains an appropriate 
model (see Fig. 1). Here M denotes the mass of a spherical shell of radius 
7' = Rq. Once again, the classical Schwarzschild limit is recovered as (r — 
Ro)/y/a — > oo, while new physics can only be expected if r — i?o ~ \fo>- 

Using solutions ([5} and ([S]) in the Einstein field equations, one can readily 
get the following expressions for p r and p t : 



1 I 



AM fZ (r-R ) 2 



Rr 1 \2' 



4q 
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(since v' = l/r) and 
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Moreover, since A' = - — r „ m 

' r r—2m 

arrive at 

Pt = 



and 
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17) 



As a final comment, the spacetime is not asymptotically flat. If this space- 
time were to be joined to an exterior Schwarzschild spacetime at the bound- 
ary of the halo, then the pressure anisotropy would be an advantage. The 
reason is that, according to Ref. [23], the solution cannot be matched to the 
Schwarzschild exterior metric at the boundary if the pressure is isotropic. 



4 Stability of circular orbits 



Given the four-velocity U a — of a test particle moving solely in the 



subspace of the halo and restricting ourselves to 9 



-mfy can be cast in a Newtonian form 



r-\ =E Z 



V(r) 



7r/2, the equation 
(18) 



which gives 



V(r) 



E 2 1 



B {] 



1 + ^r 



r 



The constants 



U 



E = — and L = — — 



m 



m 



(19) 



(20) 



are, respectively, the conserved relativistic energy and angular momentum 
per unit rest mass of the test particle [24] . Circular orbits are defined by 
r = Ro, a constant, so that = and, additionally, 4p \ r= R ~ 0. From 
these two conditions follow the conserved parameters: 



L = ± 



and, using L in V(R) = —E 2 , we get 
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(22) 



r =Ra< and unstable if ^X- \ r =R > 0. We 



The orbits will be stable if 
first determine at r = Ro by substituting the expressions for L and E: 



d 2 V 
dr 2 



= A + B + C + D, 



r=Ra 
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where 



A = 



1 - 
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(23) 



In Eqs. ([I]) and (jTDJ) , r is the distance in the radial outward direction. We 
will therefore assume that r > Rq, which does not result in a loss of gener- 
ality. Our aim is to show that under certain conditions, V" (Rq) is negative, 
resulting in a stable orbit. 

Since a controls the "width" of the Gaussian curve, we study the smearing 
effect by assuming that r — Rq ~ y/a, the fundamental condition discussed in 
Sec. [3j Returning to Eq. (|23l) . the mass of the shell, which we are treating as 
the analogue of a smeared central object, has a relatively small value when 
measured in kpc (and we must therefore assume that G = 1). In conjunction 
with the condition r — Rq ps s/a, the first term A is less than I/Rq. Next, 
factoring B, we have 



r(r-R ) 2 
Act 



(r - Rq) 



2M r 
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/7ra 3 / 2 
Rq 



R c -(r-R a ) 2 /4a 
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Since (t — Rq) 2 w a, the first factor is approximately equal to Ro — r/A, which 
is positive for r not too large. Observe that B exceeds C in absolute value. 
Because of the a in the denominator, \B\ easily overtakes A as well, even 
for moderately small a. Finally, since D is negative, V"(Rq) is also negative, 
provided, of course, that Rq — r/A is positive. The last condition limits the 
size of the variable r. 
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4.1 The effect of the noncommutative geometry 

It is shown by Nicolini, Smailagic and Spallucci [2T] that at large distances 
one would expect only a minimal deviation from the standard vacuum Schwarzschild 
geometry. Referring to Eqs. ((J) and (fTU)) . new physics can only be expected 
if r w \fa. In our model, this condition corresponds to r — Rq ps </a. As we 
just saw, if r is too large, then V"(Rq) is no longer negative. This case is an 
example of a large distance that produces the appearance of a Schwarzschild 
geometry since the smearing, although very much present, is no longer ap- 
parent. The reason is that, given the present state of research in this field, 
it is not clear whether the value of alpha can be determined. An interesting 
consequence of this lack is that, in a sense, the unseen dark matter is replaced 
by the unseen noncommutative geometry. (Regarding large-scale structures, 
it is worth noting that a Gaussian source has also been used in Ref. [27] 
to model phantom-energy supported wormholes, as well as in Ref. [2H] to 
model the physical effects of short-distance fluctuations of noncommutative 
coordinates in the study of black holes). 

Let us now come to the issue of reproduction of flat rotation curves from 
our model. It can be observed that the equation §5§ has also been employed 
by Farook et al. [TS] (vide the equation (9) therein). Therefore, based on the 
discussion made in the beginning of the Sec. 3 of the present investigation and 
following Farook et al. |15) , we can reproduce the behaviours of the present 
model in the Fig. 2. It is just a graph to illustrate that the correct shape can 
be achieved by using the following ansatz v§ of the form = ar exp(— k\r) + 
/3 [1 — exp(— k2r)\. Note that the constants a and j3 in this ansatz can have 
any units we want, and hence the whole graph can be arbitrarily re-scaled. 
A simpler way to put it is that the parameters a and j3 in the ansatz can 
be chosen to have values that fit the observational data. If one wants the 
graph to be quantitatively correct (not just qualitatively) one just needs to 
let a — > a/3 and (3 — > (3/3, but we did not use the data for Fig. 2 in the rest 
of the paper. Just re-scaling the parameters, one gets the desire figure. In a 
similar way the behaviours, via Fig. 3, can be seen for e u vs r. 



4.2 Negative radial pressure 

It is emphasized in Ref. |21) that a negative radial pressure is needed to 
retain the smearing effect near the origin, referring now to Eqs. ([T]) and 
(1101) . The negative pressure counteracts the inward gravitational pull, thereby 
preventing a collapse to a matter point. 

Returning to Eq. (fl5|) . if r — Rq sa \fa, then p r is indeed negative if 
M, the mass of the shell, is sufficiently large (since the shell need not be 
arbitrarily thin). But as long as r — Rq ~ y/a, we have V"(Rq) < 0, so that 
the negative pressure is attributable to the noncommutative geometry, much 
like traversable wormholes sustained in this manner |25) . We will return to 
this observation in Sec. [6l 



0.002 



velocity 



0.001 




Fig. 2 A sample rotation curve utilizing an ansatz equation for vj, of the form 
= ar cxp(— kir) +/?[! — exp(— J(2r)]. 




1.000000 ^ — ' ■ ' ■ ' ■ 1 

5000 10000 15000 20000 

r(pc) 

Fig. 3 Plot for the variation of e" vs r. The upper and lower panels correspond to 
short and long r behaviours. The dotted, dashed and long-dashed curves represent 
v$ — 200, 250 and 300 Kms/second, respectively. For all these, thin curves and 
thick curves represent Bo — 0.9999995 and Bo = 1.0, respectively. The chain, solid 
and thick solid curves respectively, represent = 200, 250 and 300 Kms/second 
but for Bo = 1.0000005. 
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5 Attraction and total gravitational energy 

Our next step is to consider the question of attractive gravity by studying 
the geodesic equation for a test particle that is moving along a circular path 
of radius r = R: 



dr 2 ^ dr dr 



This equation implies that 



(24) 



dr 2 







(dt\ 






\&t) 



(25) 



making use of Eq. ©. As before, as long as ^- — 0, we get 



a r 
dr 2 



1 x ^. , ; i ( dt \ 



dr <°- 



(26) 



We conclude that objects are attracted toward the center. 

According to Lyndell-Bell et al. [2(5], we can determine the total gravita- 
tional energy Eq between two fixed radii r\ and by means of the following 
formula: 



E G = M N - E M = 4tt f [1 - Ve x ( r )]pr 2 dr, 

Jri 



(27) 



where 



Mat = 47r / pr 2 dr 



(28) 



is the Newtonian mass given by 



Mat = 47r / pr dr 
J n 
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(4 7 r) 1 /2( Q ,)i 
So the total gravitational energy is 
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(4tt) 1 /2( c 



rr 2 


1 - 


_^ 2m(r) 
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(r - R f 
4a 



r 2 dr. (29) 



(r - Rpf 
4a 



r 2 dr, 



(30) 

where m(r) is given in Eq. (fTB")) . As noted in Sec. 4, M is relatively small. So 
the integrand, and hence Eg, are negative, showing that gravity in the halo 
is indeed attractive. 
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6 The observed equation of state 

The equation of state of the halo fluid can be obtained from a combination 
of rotation curves and lensing measurements. To this end, let us rewrite the 
metric, Eq. in the following form: 

ds 2 = _ e 2*(r) dt 2 + 1 -—dr 2 + r 2 + r 2 (d9 2 + sinHdtf), (31) 

2m(r) j v 7-/7 \ / 

where 

<Z>(r) = - [In B + I In r] (32) 

and m(r) is given in Eq. (TTBl . As discussed in Ref. [23], the functions are 
determined indirectly from certain lensing measurements defined by 

<2>(r) 1 fm(r) lnB \nr l f M 7 (3/2, (r — R ) 2 /Aa) 

—^dr — — ■ 1 ■ h / =-= -dr 



2 



and 



(33) 



1 ixj, s 1 ^ M 7 (3/2,(r-i? ) 2 /4a) 
m, enfl = -r 2 <?'(r) + -m(r) = - + V ^- (34) 

The observed equation of state depends on the dimensionless quantity 

3 P 3 2m ; e „ s - m RC 

due to Faber and Visser [3U]. The subscript i?C refers to the rotation curve, 
i. e. 

<l>RC = ${r) = ^[lnB + nnr] (36) 

and 

turc = r 2 <P'(r) = |. (37) 
The prime denotes the derivative with respect to r. The result is 

_ 2 m' RC - m' lens _ U/nai - M(r - R ) 2 exp [-(r - Ro) 2 /4a] 



3 Mens ~ m RC 3M ( r " ^o) 2 exp [-(r - i?o) 2 /4a] 

(38) 

At first glance, tu(r) > due to the small value of M, just as in the 
case of ordinary matter. But if a is very small and M sufficiently large, a 
situation we encountered in Sec.|4j then u)(r) < 0. This also follows from the 
fact that pt < for a sufficiently small a [Eq. (|17|)]. since we already know 
that p r is negative. So the "quintessence-like" condition is due entirely to the 
noncommutative geometry. Observe, however, that as a gets very small, uj(r) 
approached —1/3 from the right, so that w(r) is not in the actual quintessence 
range. 

As a final comment, at large distances, where the pressures are no longer 
negative, we have w(r) > 0. This may be interpreted to mean that the smear- 
ing, although still present, is no longer seen at this distance, resulting in the 
appearance of ordinary matter. 
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7 Conclusion 



In this paper we assume the existence of fiat rotation curves in galaxies and 
discuss such curves in the context of noncommutative geometry. The energy 
density of this geometry is a smeared gravitational source, where the mass 
is diffused throughout a region of linear dimension \fa. This diffusion is 
due to the uncertainty attributable to the commutator and is applicable to 
other large-scale structures. It is shown that this intrinsic property of the 
noncommutative geometry is able to account for the stable circular orbits in 
remote regions of the halo, as well as for attractive gravity. Sufficiently far 
from the orbit, the smearing is no longer observed, which may be interpreted 
to mean that stable orbits are due to the unseen noncommutative geometry 
instead of the unseen dark matter. 
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